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ABSTRACT 

High-speed spectroscopy of two pulsating subdwarf B stars, KPD 2109+4401 and 
PB8783, is presented. Radial motions are detected with the same frequencies as re- 
ported from photometric observations and with amplitudes of ~ 2kms~ 1 in two or 
more independent modes. These represent the first direct observations of surface mo- 
tion due to multimode non-radial oscillations in subdwarf B stars. In the case of the 
sdB+F binary PB 8783, the velocities of both components are resolved; high-frequency 
oscillations are found only in the sdB star and not the F star. There also appears to 
be evidence for mutual motion of the binary components. If confirmed, it implies that 
the F-type companion is > 1.2 times more massive than the sdB star, while the am- 
plitude of the F star acceleration over 4 hours would constrain the orbital period to 
lie between 0.5 and 3. 2d. 

Key words: stars: individual: KPD 2109+4401, stars: individual: PB 8783, stars: 
oscillations, stars: variables: other (EC14026/sdBV) 



1 INTRODUCTION 

The discovery of small-amplitude non-radial pulsations in a 
number of subdwarf B stars (sdBVs) has introduced a pow- 
erful new tool for the study of stellar remnants (Kilkenny et 
al. 1997=SDBVI and subsequent papers). Long time-series 
photometric campaigns have detected rich spectra of oscil- 
lations in over a dozen targets, with frequencies and ampli- 
tudes generally indicative of low-order (1 = — 2) and low- 
degree modes. The discovery has revolutionized the study 
of subdwarf B stars for the simple reason that pulsations 
in these stars were not expected. From a theoretical point 
of view, it appears that the pulsation mechanism is only ef- 
fective when diffusion processes modify the outer layers of 
the star. Metal-enrichment in a specific layer must conspire 
with the local temperature to drive pulsations through the 
opacity mechanism (Charpinet et al. 1997). 

In terms of their effective temperature and surface 
gravity, sdBVs are indistinguishable from their non-variable 
counterparts - there is no instability strip in which all sdBs 
pulsate. The origin of all sdBs remains a puzzle, although 
the increasing detection of sdB binarity, including several 
sdBVs, may point to a previous phase of common-envelope 

* Based on observations obtained with the William Herschel Tele- 
scope, ING, La Palma, Spain 



evolution in many, if not all. The structure of sdBs is par- 
tially hidden as a consequence of atmospheric diffusion, 
which disguises the true composition and the mass of the 
hydrogen-rich envelope, both of which are key diagnostics of 
previous evolution. By enabling an exploration of the com- 
position and mass of these outer layers, asteroseismology 
may be of a pivotal importance. 

Immediately after the discovery of sdBV pulsations, 
we recognized the potential for spectroscopy to provide ad- 
ditional diagnostics. Mode identification from photometry 
alone is challenging, and could be assisted by the identifica- 
tion of line-profile variations, whilst the comparison of ra- 
dial and light amplitudes could be used to determine stellar 
radii directly. Spectroscopy might also demonstrate the pres- 
ence of higher-order modes not detected photometrically. 
The limitations are that the periods are short (100-500s) 
compared with conventional CCD readout times, the stars 
are faint (12-15 mag.), and the photometric amplitudes are 
typically only a few tenths of one per cent. High-resolution 
high-S/N multi-line studies such as those obtained for non- 
radial oscillations in rapidly rotating bright O and B stars 
(e.g. Reid et al. 1993, Telting, Aerts & Mathias 1997) would 
not appear to be feasible. However, the development of new 
techniques offered the possibility to acquire high-speed spec- 
troscopy of sdBVs and here we report our first success- 
ful observations obtained in 1998 October. Subsequently, 
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Table 1. Summary of 1998 October WHT driftmode observations 
of sdBVs 
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O 'Toole et al. (2000) announced preliminary results of inde- 
pendent radial velocity observations. They report a 9kms _1 
amplitude periodic variation in the large amplitude sdBV 
PG1605+072 at the principal frequency of 2.10 mHz found 
photometrically by Koen et al. (1998, SDBV VII). 



2 OBSERVATIONS 

The targets selected for our initial study were PB 8783 
(Koen et al. 1997 = SDBV II) and KPD 2109+4401 (Koen 
et al. 1998 = SDBV XI, Billeres et al. 1998). Although not 
ideal, our observations were scheduled when possibly more 
exciting targets such as PG1336-018 (Kilkenny et al. 1998 
= SDBV VIII) and PG1605+072 (SDBV VII) were not ac- 
cessible. On the other hand, both targets have been well- 
studied photometrically (O'Donoghue et al 1998 = SDBV V, 
SDBV XI, Billeres et al. 1998) and show multiple periods in 
the ranges 7-10 and 5-6 mHz respectively, with amplitudes 
of a few millimagnitudes. 

Observations were obtained with the 4.2m William Her- 
schel Telescope on 1998 October 3 and 4 using the blue 
arm of the intermediate dispersion spectrograph ISIS. The 
R1200B grating and the TEK1 CCD with 1124 2 24^m pix- 
els yielded an instrumental resolution (2 pixel) R = 5 000 
in the wavelength interval 4020-4420 A. This wavelength 
region was chosen because it contains two strong Balmer 
lines and, potentially, a number of neutral helium and mi- 
nor species lines that are normally observed in early-B stars. 
It also maximises the photon collection rate. 

The CCD was read out in low-smear drift mode (Rutten 
et al. 1997), in which only a small number j of CCD rows 
(parallel to the dispersion direction) are read out at one 
time. A dekker is used to limit the slit-length, thus only a 
fraction of the CCD window, ~ j rows, is exposed at one 
time. After exposing for an interval t, the CCD contents 
are stepped down by j rows. This charge transfer imposes a 
small dead time d which is roughly proportional to the CCD 
window size. Each set of j rows is accumulated into a data 
cube containing n individual 2D spectra, stacked adjacent 
to one another. In practice the first few spectra in the frame 
are null, corresponding to the unexposed part of the CCD 
being read out before exposed rows reach the CCD edge. 
Comparison lamps and stellar spectra can be obtained in 
exactly the same way. Maximum frame sizes limit n to a few 
hundred; we normally adopted n = 200 before beginning a 
new frame. 

In our experiment, j was set just wide enough to admit 
a small region of sky either side of the stellar image and 
to allow for edge effects between adjacent spectra. On-chip 
binning by a factor of 3 further reduced the read-out time 
and data-volume, t was chosen so that t + d was approxi- 
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Figure 1. The mean spectrum obtained for each of two pulsating 
sdB stars observed in 1998 October. See text for further details. 



mately 1/10 of the typical observed pulsation period. Values 
for t, j and En are given for each target in Table [j], together 
with the average number of counts detected per wavelength 
resolution element c, and the maximum S/N ratio s antic- 
ipated in each exposure from nominal performance figures 
(SIGNAL, Benn 1997). In this configuration, d = 1.57s and 
the read-out noise was 5.24 electrons per pixel. 

Due to occasional errors in the CCD controller, not all 
exposure times t were equal, roughly one in 100 exposures 
would drop by an unpredictable amount. These were identi- 
fied by integrating the total counts in each individual spec- 
trum and adjusting individual exposure times to reflect the 
photometry. 

The data cubes were reduced using scripted standard 
IRAF routines. After bias subtraction and natfielding, the 
1-dimensional spectra were optimally extracted and sky sub- 
tracted. The copper-argon arc spectra obtained at the be- 
ginning and end of each data cube were time weighted for 
each stellar observation and the corrected calibration de- 
rived and applied. Over the duty cycle of one data cube the 
arc shifts were never greater than a few microns. Further re- 
ductions included normalization of the spectra and cosmic 
ray removal. 

Cross-correlation of the arc spectra showed that the 
maximum shift during a half hour period was around 7 mi- 
crons (about a quarter of a pixel). Arc exposures showed 
that this was varying smoothly during the night and as cali- 
bration data were obtained before and after each stellar data 
cube, then simple spectral interpolation of the extracted one 
dimensional spectra (weighted by the time interval from the 
bracketing arc spectra) removed any instrumental shift to 
probably within 1 micron (< 0.9kms _1 ). 

The final data product is a set of files each contain- 
ing n 1-dimensional wavelength-calibrated flux-normalized 
spectra. Each spectrum is time-tagged. The mean spectrum 
for each of our targets is shown in Fig. |l| 



3 ANALYSIS AND RESULTS 

The primary objective of our observations was to detect 
and measure radial motion due to the dominant oscilla- 
tory modes detected from photometry. There was the ad- 



© 0000 RAS, MNRAS 000,|l|-| 



Radial velocities of pulsating subdwarf B stars 3 



KPD 2109+4401 





60 




40 


\ 




a 


20 


^_ 




>. 









;j 
o 


-20 






> 


-40 




-60 



.35 



.40 .45 

JD - 2451091 



.50 




10 20 30 40 50 60 70 80 90 
Frequency (mHz) 




4 5 6 

Frequency (mHz) 

Figure 2. Radial velocities and amplitude spectrum of 
KPD 2109+4401 from WHT high-speed spectroscopy on 1998 Oct 
4. The slowly varying function b(t) is shown in the top panel. The 
bottom panel is an enlargement of the marked section of the am- 
plitude spectrum shown in the middle panel. The amplitudes of 
photometric frequencies inset in the bottom panel are are taken 
from SDBVIX. 



ditional possibility that line-profile variations or multiple 
modes might be detected, as well as any mutual motion 
within any binary system. 

The most commonly adopted method for measuring 
radial-velocity shifts in stellar spectra is the cross-correlation 
function (ccf). We constructed a mean spectrum for each 
target. This template was cleaned to remove stationary fea- 
tures such as bad CCD columns. Both template and in- 
dividual spectra were continuum-subtracted; they were al- 
ready linearized to the same wavelength scales. The ccf for 
each spectrum relative to the template was computed and 
stored as a 2d function v), where t is now time, St is 
the mean sampling interval and v is radial velocity relative 
to the template. The velocity of each spectrum was found 
by a) fitting a gaussian to, b) fitting a parabola to and c) 
computing the centroid of x( v ) m t ne range — Sv to +8v, to 
give functions v g (t),v p (t), and v c (t), respectively. Moments 
Mj(t) — J(v — v TC fyx{ v ,~t)dv;j — 1,2,3 were also com- 
puted. With ti ro f = 0, Mi corresponds to the centroid. v(t) 
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Figure 3. As Fig. | for PB8783 (1998 Oct 3). The photomet- 
ric frequencies are taken from SDBVV, mmi denotes parts-per- 
thousand and 1 mmi equals 1.086 mmag. 



was sensitive both to the method chosen and the value of 
Sv. 

The functions v(t) were investigated for periodic con- 
tent by means of the discrete Fourier transform T . Either 
a low-degree polynomial or a long-period sine function b(i) 
was first fit and subtracted from v(t). J-{v(t) — b(t)} was 
then computed. 

The velocity data v p (t) and T are shown in Figs. ^| and 
^| These figures also show a representation of the princi- 
pal frequencies and amplitudes discovered photometrically 

V (SDBVV,IX)Q 

The information we require from the amplitude spec- 
trum is not in the first instance the significance of the peaks 
but rather the amplitude associated with expected frequen- 
cies. Over the entire range 0.1 ~ 2fSnSt < v/mHz < 1/St ~ 
90, T shows many peaks of comparable amplitude which on 
their own have little statistical significance. However the co- 
incidence in Figs. ^ and ^ of the highest peaks in T with the 
highest peaks in V is more remarkable and is regarded as 

t Throughout this paper, the term amplitude applied to peri- 
odic signals refers to the semi-amplitude a of the sine function 
a sm(2iTut + <f>) 
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Table 2. Velocity amplitudes of detected pulsation modes in 
KPS 2109+4401. Up is the photometric frequency, a p the photo- 
metric amplitude (SDBVIX), v v the frequency of the measured 
velocity peak, a v the measured velocity amplitude. 
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Table 3. Simulated multifrequency observations: u, a are the fre- 
quencies and amplitudes of the test signal, n is the number of test 
runs, a is the standard deviation of normally-distributed noise, a' 
are the mean and variance of the highest detected signal within 
±0.1 mHz of the principal test frequencies. 
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1.30 
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1.96 ± 0.30 


0.92 + 0.15 


7.862 


1.40 
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0.82 + 0.13 
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2.40 
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10.619 
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a positive signature of oscillatory behaviour in the spectro- 
scopic data. The resolution 2/n8t ~ O.lmHz in T is insuf- 
ficient to resolve the fine structure observed in much longer 
photometric time series. 



3.1 KPD 2109+4401 

Because its spectrum is not contaminated by that of a late- 
type companion, the results for KPD 2109+4401 are the less 
complicated to interpret. A background drift of approxi- 
mately 9 kms -1 over 5 hours is approximately linear and 
defines b. We adopted Sv — 150 kms -1 to establish v p (t) 



and to construct Fig. |2| Principal peaks in the amplitude 
spectrum are evident at v\ = 5.09 mHz and V2 = 5.48 mHz. 

The same frequency structure in T was obtained for 
Sv = 150,300 and 450 kms' 1 when either a gaussian, 
parabolic fit or centroid is used to measure v(t). In the in- 
terval 4.5-6.4 mHz, minor peaks occured at 1^3-7 = 4.77, 
5.24, 5.33, 5.90 and 6.03 mHz in most of these analyses. 
Beating may be responsible for some peaks; for example 
f3 = v\ — (u2 — vi), ve = V2 + (y2 — fi). Thus only v\ and 
V2 could be identified as genuine from the present data. 

The probability of finding a set of peaks at prespecified 
frequencies and exceeding a given amplitude threshold by 
chance may be computed as follows. For a given frequency 
range Av, count the number of peaks n which exceed a given 
threshold Vthresh- Then count the number of peaks n-^ which 
lie within a resolution width ±Sv of the predicted frequen- 
cies. The probability of this number occurring by chance is 

25v 



Av 

For KPD 2109+4401, we chose Av = 6 mHz, v thlash = 
1.5 kms -1 and Sv — 0.05 mHz which gave n — 3, n\d = 2 
and hence p — 0.0025. 

In order to interpret the amplitudes of the identified 
peaks, we carried out the following numerical experiment 
(table ^). We first formed the sum of a series of sine functions 
each having a frequency Vi, amplitude ai and a random phase 
shift (pi. The i/j and ai were based on the photometric results 
for KPD 2109+4401 (SDBVIX). This sum was sampled at 
the same times tj as the observations, 

s(tj) = Eidi sin(27r^i£j + </>i). 

To this we added normally-distributed noise with a stan- 
dard deviation a — 15kms _1 similar to the observa- 
tional scatter. T was derived and usually found to con- 
tain peaks at the principal input frequencies, but not al- 
ways with the same amplitude ratios. Minor peaks were less 
frequently identified. The amplitudes of the highest peaks 
within ±0.1mHz of the test frequencies were measured. This 
procedure was repeated n times and the average amplitudes 
for each detected frequency {a[) were computed. In the case 
of KPD 2109+4401, the spacing of some input frequencies 
is less than the resolution provided by the limited time se- 
ries. Consequently, the measured amplitudes represent some 
sum from several independent sine functions. Interference, or 
beating, means that the apparent amplitude of such a signal 
will vary according to the relative phases of the components. 
This is apparent in table ^ (six- frequency simulation), where 
the relative standard deviation of the amplitude of the com- 
posite signal at 5.09 mHz is much higher than that of the 
single- frequency signal at 5.21 mHz. The experiment was re- 
peated with only two well-resolved input frequencies and 
yielded the ratio a' /a ~ 1.0 ± 0.1. Reducing a reduces the 
error on a' /a due to experimental noise, but not that due 
to mode beating or to data sampling. It was also repeated 
for v and a appropriate to the second target PB 8783 with 
a similar result. 

Therefore, we conclude that the peak amplitudes of 
2.49 and 2.68 kms -1 at frequencies 5.09 and 5.48 mHz in 
KPD 2109+4401 reflect to within ~ 20% (standard devi- 
ation) the combined real amplitudes of oscillations within 
~ 0.05mHz of these frequencies (Table ^|). The principal 
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Figure 4. Peak of a typical cross-correlation function x(v) f° r 
PB 8783 (crosses). The gaussian + quadratic fit used to construct 
the two-component velocity solution (Figs. [] and ^) is also shown 
(solid line). 



causes of uncertainty in the amplitudes are beating between 
unresolved modes and experimental noise. 

The 9 kms -1 drift in b may be real. Instrumental drift is 
limited to < 0.9 kms -1 (section 2). The heliocentric correc- 
tion, which was not applied to the data, changes by less than 
0.4kms -1 . In the simulations described above, the r.m.s. 
amplitude of b was 1.7kms _1 . Although KPD 2109+4401 is 
not known to be a binary, it cannot be ruled out. PG 1336- 
018 was only recognised as a binary from its eclipses; 
the secondary is too faint to be detected spectroscopically 
(SDBVVIII). 



3.2 Line profile variations? 

In addition to radial-velocity information, \ contains infor- 
mation about line-profile variations, since the ccf represents 
an average line profile within the region of spectrum ana- 
lyzed. In the present case, the average line profile is dom- 
inated by two Balmer lines, H7 and H5. The average line 
width and asymmetry are indicated either by the moments 
A^2,3, respectively, or by the coefficients 52,4 in the gaussian 
+ linear fit g(v) = g exp((v - g\f I gi) + gs + g±v. 

Neither the Fourier transform J r {M2(t)} or J-{g2 (t)} 
showed peaks near v\ or V2- With suitably chosen Sv, both 
J-{Ms(t)} and J-{gi(t)} showed frequency structure similar 
to J-{v(t)}. However, the amplitude of peaks in jT{Mi(t)} 
only approached that of ^F{v p (t)} if Sv was sufficiently large, 
in which case the periodic content of M^(t) was essentially 
lost. This implied that M3 contained velocity information 
not properly incorporated into Mi. 

jT{g 4 (t)} showed peaks at 5.09 and 5.55 mHz with 
Sv ~ 400 — 600 kms -1 , but not otherwise. In the case of 
intrinsically sharp lines, such asymmetry may reflect differ- 
ential radial motion across the rotationally broadened line 
profile. In the present study, any rotational broadening is 
convolved with Stark broadening in the Balmer lines and is 
difficult to interpret. Consequently, whilst we suspect that 
there may be evidence for line-profile variations with simi- 
lar frequencies to the radial motions, we can draw no firm 
conclusions from the present data. 



Table 4. Velocity amplitudes of detected pulsation modes in 
PB8783. Symbols are as in table ^ except that v p ,a p are from 
SDBVV. The identification marked (:) is tentative. 
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3.3 PB8783 

The binary PB 8783 is altogether more difficult to analyze 
because of the superposition of sdB and F star spectra. The 
results shown in Fig. are based on fitting parabola to x( v ) 
with Sv = 420 kms - . In deriving the amplitude spectrum 
in Fig. ^| the background function b(t) was found to be best 
given by a sinusoid with an amplitude of 10 kms -1 , rep- 
resenting a systemic drift substantially larger than that in 
KPD 2109+4401. 

On inspection, the typical ccf for PB 8783 consists of a 
strong broad component (FWHM ~ 1400 kms -1 ) superim- 
posed by a small narrow component (FWHM ~ 200 kms -1 , 
Fig. Q). The narrow feature disappears when the template 
is replaced by a theoretical sdB star spectrum containing 
H and He lines only, and is deduced to be due principally 
to metal-lines in the F-type secondary. This is confirmed 
when x i s calculated using only wavelengths between H5 
and H7; the broad component disappears leaving a small 
narrow peak. 

The detection of two components in x offers the pos- 
sibility to resolve the radial motions of the binary compo- 
nents. The mutual motion of the F-star is clearly apparent 
when x excludes the Balmer lines, and appears as a change 
of ~ 10 kms -1 during the five- hour observing run. However, 
the sign of the change remains the same when the theoreti- 
cal sdB spectrum is used as a template, indicating that the 
Balmer lines in the F-star contribute significantly to x- 

The problem was resolved by fitting x{~ $ v '■ +Sv),8v = 
280 kms -1 with the function g(v) = goexp((v — g\) 2 j g^) + 
g-i + giV + g^v 2 , where the quadratic terms fit the broad com- 
ponent and the gaussian terms represent the narrow compo- 
nent (Fig. ^|). The functions gi(v) and — g^{v)/2gs(v) were 
examined to determine the mutual motion of each star and 
to look for high-frequency content. The results are shown 
in Figs. ^ and ^ where it will be seen that the slowly vary- 
ing components are of opposite sign. It is suggested that, in 
the absence of viable alternatives for these systematic drifts, 
these slowly varying components are most probably due to 
the mutual motion of the two stars within the binary. 

If this is the linear regression on the first 4 hours 

of both sets of velocity data gives the relative acceleration of 
the two components. The result, Awp/Au s dB = —0.81 + 0.10, 
provides an estimate for the mass ratio of the two stars 
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Figure 5. Radial velocities and amplitude spectrum of the sdB 
star in PB 8783. The least squares solution to the slowly varying 
velocity component is shown in the top panel (solid line). The 
amplitudes of photometric frequencies are as in Fig. til 



Mgds/Mp ~ 0.8. This is a little high for an F-star on 
the main-sequence (1.1 — 1.8 Mq) and a canonical sdB star 
(~ 0.5 Mq), and suggests that the relative motions have 
still not been fully resolved. Further information concerning 
the velocities of the two components relative to one another 
is convolved into the composite Balmer profiles in the ccf 
template and will require observations of the complete bi- 
nary orbit to disentangle. Whilst the binary period is clearly 
longer than the data sequence, it is constrained by the rela- 
tion 



K F = 210 



1 /M\ J 
1 + q \~Pj 



sin i kms 



(1) 



M is the total mass of the system (Mq), q = Mf/M^b 
is the mass ratio, P is the orbital period (d) and Kf is the 
velocity semi-amplitude of the F star ( kms -1 ). For example, 
with Af sdB = 0.55 M Q , M f = 1.2 M Q , i = 30° and P = 
2d, .Kf = 32 kms" 1 . The maximum change in velocity over 
the timescale of our observations (0.15d) would be Auf ~ 
15 kms" 1 , compared with 13 kms" 1 observed. Inverting the 
procedure, we can determine the maximum P as a function 
of i required to obtain Avf = 13 kms" 1 in 0.15d. For Mf = 
1.2 — 1.8 Mq, the maximum P increases from 1.0 — 0.8d at 
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Figure 6. As Fig. § for the F star in PB 8783. 



i = 10° to P = 3.7 - 3.2d at i = 80°. Eclipses are not 
expected at inclinations < 80° . 

It is satisfying to detect the photometric frequencies 
already reported for PB8783 (SDBV) in the deconvolved 
radial velocities of the sdB star (Fig. H) and to find no such 
signal in the radial velocities of the F star (Fig. |^). This 
provides confirmation that the photometric and radial ve- 
locity variations are due to pulsations in the sdB star and 
not in the cool companion. A similar result was obtained 
for PG1336-018 (SDBV VIII), except in that case the per- 
sistence of oscillations during an eclipse of the unseen com- 
panion provided the necessary proof. 

From this point, the interpretation of the sdB star ve- 
locities is the same as that for KPD2209+4101. Peaks in 
the velocity amplitude spectrum were identified at five fre- 
quencies also present in the photometry. The resulting pul- 
sation velocity amplitudes were thus found to be in the 
ranee 1.1 - 2.2kms _1 and are shown in Table ^. Regard- 
ing the probability that these frequencies were identified 
by chance and following the prescription already given for 
KPD 2209+4101, we found with Wthresh = 1.5 kms" 1 that 
n = 9,n id = 3 and p = 0.0034 (Fig. ^). Relaxing Wthresh = 
1.0 kms" 1 yielded n = 27, n id = 6 and p = 0.0083. 
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3.4 Relating luminosity and velocity amplitudes 

Kjeldsen & Bedding (1995) investigated the relation be- 
tween velocity and luminosity amplitudes for several classes 
of main-sequence, giant and supergiant pulsators and de- 
rived the calibrated relation 

^ = (A/550^(W5777IO a 20 - 1PPm| ^ 
where ppm denotes parts-per-million. 

We tested whether these predictions could be extended 
to the sdBVs since we have both light (SL/L) and velocity 
(i> osc ) amplitudes for the same oscillation frequencies. We 
found that the observed velocity amplitudes are > 2 times 
those expected from the photometry and Eqn. (S). However 
Eqn. has been derived from cool stars pulsating in radial 
or non-radial p-modes and it may not be applicable to pul- 
sations in the hot sdBVs. We note that the amplitudes of 
individual modes in sdBVs have been reported to vary from 
one season and possibly one night to another (SDBV/,VIII). 
Beating between unresolved modes will certainly affect the 
velocity amplitude measurements. Simultaneous high-speed 
photometry and spectrocopy will be required to properly re- 
late the velocity and luminosity amplitudes and to make any 
futher validation of Eqn. ^. 

3.5 For future investigation 

We have left two issues for future investigation. 

1) Compared with its apparent surface gravity, the observed 
pulsation frequencies of PB 8783 are too high (SDBV XII) . 
The expectation is that its surface gravity (logg = 5.55, 
O'Donoghue et al. 1997, SDBV IV) should be much higher. 
It is recognised that this could be due to inadequate sub- 
traction of the F-star spectrum, a conjecture which should 
be checked using higher quality data, including those pre- 
sented here. 

2) If the mode and amplitude of a spherical harmonic re- 
sponsible for each pulsation frequency is specified, then it is 
relatively straightforward to integrate the projected surface 
velocities to provide a theoretical radial velocity curve. It is 
slightly more complicated to construct the theoretical light 
curve. From the combination it should be possible to refine 
the mode and amplitude measurements already established. 



4 CONCLUSIONS 

High-speed high-resolution spectroscopy of two non-radially 
pulsating subdwarf B stars KPD 2109+4401 and PB 8783 
has been acquired with the William Herschel telescope. 
These data show radial velocity variations at both high- 
and low-frequencies. High-frequency peaks in the velocity 
amplitude spectrum correspond to frequencies identified in 
the light curves of both stars, and have allowed the velocity 
amplitudes associated with these stellar oscillations to be 
estimated. Typically ~ 2kms _1 , these translate into radial 
variations of some 100 — 250km within 60-90s. The prospect 
for measuring the velocity amplitudes of non-radial pulsa- 
tions in sdBVs with much higher amplitudes is promising. 
Line profile variations were not detected in the data al- 
though there are good reasons to suppose they should be 
present. The possibility that frequencies not present in the 



photometry may contribute to the velocity amplitude spec- 
trum should also be pursued. 

In the case of the binary sdB PB 8783, low-frequency 
velocity variations corresponding to the mutual motion of 
both components have been resolved. These provide an up- 
per limit to the orbital period of between 0.8 and 3d, de- 
pending principally on the orbital inclination. Establishing 
the orbital period for this star should therefore be a prior- 
ity for future observations. The deconvolution of component 
star velocities also demonstrated that high-frequency veloc- 
ity variations occur only in the sdB star and not in the F 
star. 
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